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The large self-consistency programme is reviewed. As an application the QCD /3-function is 
computed at 0{1/Nf) and the anomalous dimensions of polarized twist-2 singlet operators are 
determined at the same order. 



1 INTRODUCTION 

The principle of computing higher order perturbative corrections to quantities in a quantum 
field theory is simple to state. Specifically one wishes to determine the coefficients of a series 
expansion in powers of the coupling constant, g, which is assumed to be small. The difficulty 
is that each coefficient arises from a set of Feynman diagrams. Since the number of graphs 
increases rapidly with each order in g, this becomes a huge computational exercise. Indeed 
few calculations are available beyond the third order for four dimensional gauge theories, [1- 
3]. Such computations, however, can be complemented by alternative methods. For example, 
if a quantity whose perturbative structure is required, depends on parameters other than the 
coupling constant, then it may be possible to perform some sort of expansion in those other 
variables. In other words, to approach the calculation of the coefficients from another angle. Such 
a possibility exists for models where the basic fields lie in a multiplet of an internal symmetry 
group like SU{Nf). Then an expansion in powers of the dimensionless variable 1/Nf can be 
performed when Nf is large. Included in this class of theories are QED and QCD. 



2 BASIC FORMALISM 

We can express this more concretely for such theories by considering the structure of a typical 
quantity, F{g), whose expansion is of the form. 

Fig) = aig + (aaiV/ + 61)5' + (asiV/ + feaiV/ + ci)^^ + Oig^) (1) 

The aim is to determine {oi}, {6j} and {q} and so on, with respect to some renormalization 
scheme. For this article we will consider only the MS scheme. We have included the dependence 
on Nf explicitly. Perturbative calculations systematically extract {ai}, {02,61} and {a3,b2,ci} 
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and so on. However, F{g) can equally be calculated as an expansion in powers of with 
large provided ^A^^ is held fixed. Such an approach is in effect a reordering of the set of 
Feynman diagrams making up F[g). Clearly the coefficients which are computed successively 
at each order in l/iVf are {aj}, {6,} and {q} and so on. The question that now arises is how to 
determine an infinite set of coefficients in a compact way. 

For the first part of this article we review the method to compute the structure of the 
functions of the renormalization group equation, (RGE). This is achieved by using ideas from 
statistical physics and the critical renormalization group and which were initially applied to the 
0{N) a model in Q. Basically at a phase transition universal quantities can be measured ex- 
perimentally or calculated theoretically which characterise the physics. These critical exponents 
are fundamental and are determined from the renormalization group functions evaluated at the 
fixed point. They are functions of the space-time dimension, d = 2^, and the parameters of any 
internal symmetries of the underlying field theory, For example, the exponent r/ = r](d,Nf) 
is related to the wave function renormalization. With this relation to the RGE functions, if one 
computes critical exponents directly, in say the 1/Nf approximation, and provided one knows 
the location of the critical coupling, gc, as a function of d in powers of 1/Nf, then one can deduce 
the coefficients of the RGE functions from an e-expansion of the exponent where d = 4 — 2e. 

For QED and QCD this procedure is possible provided one considers each theory in d- 
dimensions, when 13(g) has a non-trivial zero. Specifically, 



f3{g) = id- i)g + 

has a fixed point at, [H], 



(2) 



5c = 3e/[T{R)Nf] + 0(l/iV/) (3) 

where the colour group Casimirs are defined by tr(r'^r^) = T{R)6''\ T^T" = C2{R) and f^^'^f'^^ 
= C2{G)5°'^ with Tfj the generator of the colour group and /"''^ its structure constants. So if 
F{g) is the quark wave function renormalization then 



Y. "„(3£/T(iJ))" 



.n=l 



+ 0{l/Nf 



and the coefficients {a,i} are clearly related to those of ??i(e) = F{gc) where rj = Yl'^=i ^n(e)/A^/"'- 

Exponents such as r/ are defined by considering the action as a dimensionless quantity at gc- 
With the usual QCD lagrangian 

L = iij'^l/hp^' - (G^jV4e2 (4) 

where 1 < i < Nf, = + iT^Al and G^^ = d^^A^ - d^A^ + f^'^'^A^^A^, then the quark 
and gluon propagators at criticality will have the asymptotic form, Q, in an arbitrary covariant 
gauge with parameter 6, as A;^ ^ oo, 



r h..h.^ 

(5) 



(fc2)M-/3 



A;2 



Here A and B are /c-independent amplitudes and the anomalous parts of the dimensions a and 
/3 are defined relative to the canonical pieces asQ = // — l + r//2 and /? = 1 — r/ — The 
exponent x is the quark gluon vertex anomalous dimension. Indeed each operator built out of 
the fields of (4) will have an associated critical exponent. 
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Several leading order exponents have been determined in QCD. For example, in the Landau 
gauge at leading order in 1/A^/, 

' T{R)Nf ' ' ^ 2{fi - 2)T{R)Nf ^ ' 

where r/° = (2/Li — — 2)r(2;u)/[4r^(/i)r(/i+l)r(2 — //)]. These have been deduced respectively 
by substituting the critical propagators (5) into the quark and gluon 2-point functions and the 
quark gluon 3-point function and examining the scaling behaviour of the resulting integrals in 
the critical region. The expansion of (6) in powers of e agrees exactly with the explicit 3-loop 
perturbative calculation of the corresponding renormalization group functions in the Landau 
gauge, IQ. One feature which leads to a simplification in these calculations is that one does not 
need to consider any diagrams where there are radiative corrections on an internal line. The 
reason is that one is using propagators which have non-zero anomalous dimensions. As these 
represent the effect of such corrections, including them in a calculation would in effect be a 
double counting. For higher order corrections this reduces the number of diagrams that 

need to be computed at criticality. 



3 QCD /^-FUNCTION 

We now consider the computation of the QCD /J-function as an application of the above formal- 
ism. First, we need to recall another simplification of the critical point approach. Sometimes 
more than one field theory underlies the description of a phase transition. In this case there is 
a choice of models which can be used to determine the critical exponents and such models are 
said to be in the same universality class. An example of this is the equivalence of the 0{N) a 
model and 0^ theory with an 0{N) symmetry. In three dimensions the exponents calculated in 
either determine the critical behaviour in the Heisenberg ferromagnet. However, if one model 
in the universality class has a simpler structure then an exponent calculation could be substan- 
tially reduced by computing with it. This is the case for (4). In the large Nf limit it has been 
shown, 0, that QCD is equivalent to the non-abelian Thirring model, (NATM), which has the 
lagrangian, 

L = - {A-^f/2X (7) 

where A is the coupling constant which is dimensionless in two dimensions. It is more efficient 
to use (7) to deduce universal critical exponents due to the absence of the additional triple and 
quartic gluon self-interactions. It was also shown in |^, however, that these interactions are 
correctly recovered in integrating quarks out of the gluon 3- and 4-point functions. Of course, 

(3) must be used to determine the perturbative coefficients of quantities in QCD, Q. Also, as we 
are considering a non-abelian gauge theory in covariant gauges the ghost sector must be added 
to both lagrangians. At leading order in l/Nj, though, we do not need to consider contributions 
from graphs with ghosts as they are supressed by an additional factor of 1/Nf. 

To find (3(g) at 0{1/Nf) we use a scaling law between various anomalous dimensions to 
determine the non-zero exponent lo = — P'(gc)/2. It is deduced from the gluon kinetic term of 

(4) , as 

UJ = V + X + XG (8) 

where the first two terms arise from the gluon fields in the field strength composite operator G 
= (G^i,)^ which has anomalous dimension XG- It is determined by examining a Green's function 
with G as an insertion and, using (5), calculating the leading order set of contributing 

graphs. These are illustrated in figure 1. In [P, however, the structure of the QED /?- function at 
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Figure 1: Graphs for uJi. 

0{1/Nf ) was determined by explicitly carrying out the leading order bubble sum and using MS 
renormalization. This was later extended to the critical formalism in [jlO|, which corresponds 
to the calculation of the first two graphs of figure 1. Therefore for the non-abelian extension 
one need only compute the remaining two graphs which involve the colour group factor C2{G). 
Although the critical exponent associated with the /3-function is independent of the gauge pa- 
rameter, to avoid mixing with the operator (9^yl^)^, jl^, we compute lv in the Landau gauge. 
With 6 = we find, [§, 



(2m-3)(/.-3)^ 4(2^-l)(^-l) 



(9) 



where ujq = {fi — 2). The e-expansion of (9) correctly reproduces the 3-loop (gauge independent) 
result of [1-3] and allows us to deduce new coefficients at higher orders. In the notation of (1), 

04 = - [154C2iR) + 53C2(G)]/3888 (10) 
as = C2(i2)(144C(3) + 107)/15552 + C2(G')(480C(3) - 229)/31104 (11) 

where Ci'^) is the Riemann (^-function. 

We close this section by noting some technical points in the determination of XG and anoma- 
lous dimensions of other operators. The contribution from each graph is found from the residue 



of the simple pole in the regulator of the critical theory, |12|. This (analytic) regularization is 
introduced by shifting the dimension of the gluon field by an infinitesimal amount, A. In other 
words we calculate the graphs using (5) but with the replacement /? — > /3 — A. It is important 
to appreciate that the dimension of space-time is fixed in these exponent calculations, unlike in 
explicit perturbative computations using dimensional regularization. 



4 POLARIZED OPERATORS 

We have extended this analysis to other operators whose renormalization is necessary for the 
operator product expansion used in deep inelastic scattering. For instance, one can study the 
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twist-2 flavour singlet polarized operators, |13|| , 

Og = — 5e^i"^^trG/3^D''2 _ _ _ D^'^-^G^"^ 

Of = i"5V^7^7''iL»''2 _._L,Mn^ (12) 

where 5 denotes symmetrization in the Lorentz indices and subtraction of all combinations of 
traces and n is the moment of the operator. The treatment of unpolarized flavour non-singlet and 
singlet operators at 0{1/Nf) was given in [0]. We recall that the gauge independent anomalous 
dimension of these physical operators is responsible for the way in which the associated Wilson 
coefficients vary with energy. It is important to have large information available on the renor- 
malization of (12) as the matrix of anomalous dimensions, jij(g), has recently been evaluated at 
second order, []T5| , |l^ ]. A matrix of renormalization constants arises in perturbation theory due 
to the mixing under renormalization of both operators as they have the same quantum numbers 



and canonical dimensions, |14|. 



The computation of the corresponding exponent matrix, jij{gc), is by the technique outlined 
in the previous section. We substitute the operators in Green's functions and determine the 
scaling behaviour. In the large approach there is a simplification. At criticality, due to 
the field content, the operators have different canonical dimensions and therefore do not mix. 
Alternatively this can be seen by calculating "fiji^gc) directly and observing that it is triangular 
at 0(1/A^/), ig. 

In order to relate exponent results to perturbation theory one observes that the exponent 
derived from the insertion of actually corresponds to the eigenoperator of the perturbative 
mixing matrix which is predominantly fermionic in its content. In other words working at gc one 
circumvents the need to construct operators which diagonalize the RGE governing the evolution 
of the fermionic and gluonic singlet Wilson coefficients. The result of our calculation is, for even 
n, 



(„) ^ 2C2{R)7]? 
''F-i (2^-l)(^-2)r(i?) 



(a -If Kf^-'^y 

^ (/i + n- l)(^ + n-2) 



I o / iwi/ 1 , ^ It \] At(2/i + n-5)(n + 2)5(n) 

2[fi + n — + n — 2) 



(13) 



where S{n) = r(n)r(2^)/r(2^ + n — 2) and '(/'(x) is the logarithmic derivative of the F-function. 
It agrees with the two loop results of |jl5|, |l^ when "fij{g) is diagonalized. The leading order 
1/A^ value of the dimension of Oq is equivalent to the one loop value, which follows from the 
Nf structure of the entries in jij{g). For the sake of comparison we note the exponent for the 
corresponding unpolarized operator Of = i'^Si/jj'^'^ D'^^ . . . D^"iIj is, ]14[ , 



An) 



2(m - 1) 



^ , ^ (/i-l)C2(i?)r?? 
''f'l {2^i-l){pi-2)T{R)Nf 

+ 4^[V(/u - 1 + n) - V(/i)] 

/i(2// - 1 + n)F(n - l)F(2/i) 



2ii{iJ, - 1 



i2 



{n + n-l){ii + n-2) 



(// + n - l)(/i + n - 2)F(2/i + n) 

X [(n(n - 1) + 2(^ - 1 + n)f + 2n(n - l)(/i - 2)(2^ - 3 + 2n) 
+ 4(^-2)(^-l+n)] 



(14) 



The purely four dimensional object 7^ needs to be carefully treated in the large method, 
|17|. Unlike in some approaches in dimensional regularization we can use a fully anticommuting 
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7^ in d-dimensions since our regularization is achieved by a shift in the gluon dimension, 
For closed fermion loops one performs the integral in d-dimensions before projecting with, for 
example, [17|, 

tr(7S'^7"7"7^) = ^e^.ap (15) 

As a final application we have also computed the anomalous dimension of the singlet axial 
current, Os = ip^'^j^xp. The renormalization of this operator is more involved as it is not 
conserved in the quantum theory due to the chiral anomaly, |19| ]. However, in perturbation 
theory methods have been developed to obtain the three loop anomalous dimension in the 
MS scheme. (See, for example, [|20|.) Briefly this is a two part exercise. One part involves 
carrying out a standard d-dimensional renormalization using dimensional regularization. A 
finite renormalization is also required. This is determined by ensuring the one loop character of 



the operator form of the anomaly is preserved, |2C]. A similar two part approach in 1/Nf, using 
the above rules for graphs involving 7^, yields the critical dimension of O5 as, llF 



3/^C2(fi)< .... 

^^'^ = - (/. - i)r(i?) 



The e-expansion of (16) agrees with the three loop result of pl|, 2C]. 



5 STRUCTURE FUNCTIONS 



A second area of application oil/Nf methods is in the determination of the higher order structure 
of the finite parts of Green's functions or amplitudes. Not only are such calculations important 
in obtaining new coefficients of a series, they can also be used to study summability problems. 
As the 1/A^ approach is a reordering of perturbation theory in which graphs with infinite 
chains of quark loops are treated first, one could substitute such chains into the appropriate 
amplitude and carry out the renormalization to leave a finite result. In practice this will be 
tedious and inefficient. For a better approach we recall a feature of perturbative calculations. In 
dimensionally regularized calculations with massless fields, the effect of performing a quark loop 
integral is to obtain a d-dependent result with a factor which involves the momentum raised to 
a d-dependent power. In the renormalization process this latter part is important in obtaining 
the finite piece. Motivated by this and the simple structure of the fixed point propagators, we 
can obtain 0{1/Nf) information on the finite part efficiently by using the propagators, jlT 
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We stress that this approach is not at a fixed point and is closer in spirit to perturbative calcu- 
lations. The gluon propagator has an adjusted power, 5, where 5 = 4:NfT{R)g/3. Calculations 
are performed by substituting (13) into Green's functions with undressed gluon lines with the 
integration restricted to four dimensions. The resulting 5-dependent function can be expanded 
in powers of 6 and it transpires that the coefficients are in direct correspondence with those 
obtained in perturbation theory. 

Several applications are available. For instance, the finite part of the photon propagator has 



been determined in QED, |22]. In deep inelastic scattering the 0{1/Nf) part of the longitudinal 
non-singlet structure function has been studied. The Wilson coefficients at the Lth loop is, [O], 



Qong(l'5)-^) 



d6^ 



{2 - d){l - 6)x'' 



(18) 
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where S{n, S) = T{n + 6)/[T{n)T{l + 5)/(n + 1 + 6)], and x is the Bjorken scahng variable. It 
agrees with the three loop results of p^ . 

The poles in 5 at 5 = 1 and 2 correspond to renormalons in the series. Recently, these have 
been used with the naive non-abelianization procedure, |^^, to estimate the effect of higher 
twist contributions to the full structure function, ||2^ 



6 CONCLUSIONS 



We have focussed primarily on the 0{1/Nf) corrections to the RGE functions. One important 
feature of the fixed point approach, however, is that 0{1/NJ) calculations are also possible in 
d-dimensions. For example, the anomalous dimension of the electron mass in QED has been 
determined at this order, [^] . With the simplification that occurs from the equivalence of QCD 
and NATM the extension of this and other results will be possible for the non-abelian case. 
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